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Abstract 

o : 

In this letter, we analyze the ergodic capacity of bidirectional amplify-and-forward relay selection (RS) with 

Ch ' 

imperfect channel state information (CSI), i.e., outdated CSI and imperfect channel estimation. Practically, the 

I— 5 ' 

optimal RS scheme in maximizing the ergodic capacity cannot be achieved, due to the imperfect CSI. Therefore, 
two suboptimal RS schemes are discussed and analyzed, in which the first RS scheme is based on the imperfect 
channel coefficients, and the second RS scheme is based on the predicted channel coefficients. The lower bound 
of the ergodic capacity with imperfect CSI is derived in a closed-form, which matches tightly with the simulation 
results. The results reveal that once CSI is imperfect, the ergodic capacity of bidirectional RS degrades greatly, 
whereas the RS scheme based on the predicted channel has better performance, and it approaches infinitely to 
| the optimal performance, when the prediction length is sufficiently large. 
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bidirectional relay selection, imperfect channel state information, ergodic capacity 

I. Introduction 

Recently, bidirectional relay network attracts a lot of interest, because it has better spectral efficiency than 
conventional one-way relay network when adopting the network coding technique JTj]. In addition, the relay 
selection (RS) technique has been intensively researched in the bidirectional relay network, due to its ability to 
achieve full diversity with a single relay In 0, 0, the symbol error rate (SER) of bidirectional RS 

was derived in a closed-form, which verifies that RS can achieve full diversity. The ergodic capacity analysis of 
bidirectional RS was obtained in @, 0. Furthermore, imperfect channel state information (CSI), i.e., outdated 
CSI and imperfect channel estimation, has great impact on the performance of RS, which has been fully studied 
in one-way relay network, such as the SER analysis J6], Q, and the ergodic capacity analysis (H, |9). However, 
to the best of the authors' knowledge, all the previous researches about bidirectional RS all assume the CSI is 
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perfect, and the impact of imperfect CSI, such as outdated CSI, on the performance of bidirectional RS has not 
been investigated. 

In this paper, we analyze the ergodic capacity of RS in bidirectional amplify-and-forward (AF) relay with 
imperfect CSI. The system model and the imperfect CSI model are presented in Section II. In Section III, we 
discuss the optimal RS scheme in maximizing the ergodic capacity. However, it cannot be achieved practically, 
due to the imperfect CSI. Therefore, two suboptimal RS schemes are analyzed, in which the first scheme is based 
on the imperfect channel coefficients, and the second scheme is based on the predicted channel coefficients. The 
tight lower bound of the ergodic capacity for the suboptimal RS schemes is derived in Section IV. It is noted 
that the analytical expression is provided under the generalized network structure, i.e., the channel coefficients 
follow independent but not necessarily identical complex-Gaussian fading, and the correlation coefficients of 
outdated CSI are different for different channels. In Section V, Monte-Carlo simulations verify that the analytical 
expression matches tightly with the simulated results. The results reveal that the ergodic capacity of bidirectional 
RS degrades greatly once CSI is imperfect, whereas the RS scheme based on the predicted channel can compensate 
the performance loss, and it approaches infinitely to the performance achieved by the optimal RS scheme, when 
the prediction length is sufficiently large. 

Notation: (•)*, (-) T , and (-) H represent the conjugate, the transpose, and the conjugate transpose, respectively. 
E is used for the expectation and Pr represents the probability. 



In this paper, we investigate a generalized bidirectional AF relay network with two sources Sj, j = 1,2, 
exchanging information through N relays Ri, i = 1,...,N, where each node is equipped with a single half- 
duplex antenna. The transmit powers of each source and each relay are denoted by p s and p r , respectively. The 
direct link between sources does not exist due to the shadowing effect, and the channel coefficients between Sj 
and Ri are reciprocal, denoted by hji. All the channel coefficients are independent complex-Gaussian random 
variables (RV) with zero mean and variance of <j\.., and these coefficients are constant over the duration of one 
frame. 

The whole procedure of bidirectional AF RS is divided into two parts periodically: relay selection process 
and data transmission process. During the relay selection process, e.g., the Tth frame, the central unit (CU), i.e., 
Sj, selects the best relay according to the predefined RS scheme, which will be discussed in the next section. 
During the data transmission process, e.g., the (t + r')th frame, only the selected relay is used for transmission, 
and other relays keep idle until the next relay selection process comes. Due to the time-variation of channel, the 
channel at the date transmission process hji (t + t') is quite different from that at the relay selection process 
hji (r), and their relationship is modeled as ||6] 



II. System Model 




(1) 
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where Eji is an independent identically distributed RV with hji (r); pf .. = Jq (27r/d Tr'), where Jo (•) stands 
for the zeroth order Bessel function frPZl . f^.. is the Doppler spread, and Tt' is the time delay. 

The imperfect channel estimation is also considered in this paper, and the channel hji and its estimate hji 
are related by hji = hji + eji Q, in which hji and the detection error eji follow independent zero-mean 
complex-Gaussian distributions with variances of cr? and a^., respectively, and a^.. = means no estimation 
error. 

Considering the transmission via Ri, the data transmission process of bidirectional AF relay is divided into two 
phases. During the first phase, the sources simultaneously send their respective information to Ri. The received 
signal at Ri is n = ^/pJ/i^ijSi + ■sjplhtp.i^ + where h t ji = hji (r + r'), Sj denotes the modulated symbols 
transmitted by Sj with the average power normalized, and n r i is the additive white Gaussian noise (AWGN) at Ri, 
which is zero mean and variance of cr 2 . During the second phase, Ri amplifies and forwards the received signal 
back to the sources. The signal generated by Ri satisfies ij = -y^ftrj, where = (p s lt,u +Ps7t,2i + °"n) 
is the variable-gain factor [6|, and ^ t ,ji = \ht,ji\ 2 - The received signal at Sj via Ri is yji = htjiU + n s j, where 
n s j is the AWGN at Sj. After reconstructing and canceling the self-interference, i.e., y/p s p r f3i htji ht,j%Sj [3], the 
instantaneous received signal-plus-interference-to-noise ratio (SINR) at Sj via Ri is 

A^s%ji%,ji ... 

Iji = jz rr ; (2) 

hps + Vr J %ji + i>s%,ji + C 

where ip s = p 8 ja\, ip r = Pr/^% 4>s = ips + ipsiprol^, A = ipr + 3ip s ip r a^. z + ip s ip r a% 3i , c = 2ip r ip s a^. + 
if>ri)sa% t + Vv<, + 1, and {jj} = {1, 2} or {2, 1}. 

III. Relay Selection Schemes with Outdated CSI 
The ergodic capacity of the bidirectional RS is defined as PI, 151 

C = E 7lfc , 72fe | i log 2 (1 + 7i fc ) + i log 2 (1 + 72fe)| 

where k in the index of the selected relay, 7jfc can be obtained by ©, and the pre-log factor 1/2 means that the 
transmission of one data block from one source to the other occupies two phases. 
To maximize the ergodic capacity, the optimal RS scheme is 

k = arg max min j | hi j (r + r') | 2 , \h 2 i (r + r') | 2 | (4) 

where the detailed explanation of the scheme is given in Appendix A. 

However, hji (r + r') is unknown at the relay selection process, due to the time-variation of channel (0Q). 
Accordingly, the optimal RS scheme (01) cannot be achieved with outdated CSI, thus we discuss two suboptimal 
RS schemes which can be implemented practically. 



(3) 
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The first alternative RS scheme with outdated CSI is obtained by substituting the channel hji (r + r') in (0]) 
with the outdated channel hji (r) obtained at the relay selection process J2], |@], i.e., 



k = arg max min 



ha (r) , h 2 i (t) 



(5) 



The second alternative RS scheme with outdated CSI is obtained by substituting the channel hji (r + t') in 
(01) with its predicted value h p ji (r + t'), i.e., 

|2 



fc 



arg max min 1 1 i, (r + r')| 2 , \ h p ^i (r + r') | 2 | 



(6) 



In this paper, the Wiener filter [13 1 is applied for channel prediction, which is the linear optimal prediction 
in minimizing the mean square error, thus we have h p ji{r + r'^ = w^^hj,, where hji = \hji[r) ,hji[r — 
A) , • • • , hji (t— Lj{A+ A)] T contains the current and previous {Lji — 1) channel coefficients, which are obtained 
from previous relay selection processes, A is the interval of adjacent relay selection processes in frames, Lji is the 
prediction length, and w optj j = R^ 1 ^, in which Hji = E |hjjh^| and rji = E jhjj/i*^ j. According to lfl3l . 



hji (r + r') = h Pt ji (r + r') + \Jv\.. — a Pji n ji' where h p ji (r + r') is a complex-Gaussian RV with zero mean 
and variance of a p . z = r^Hj^Tji, riji is an independent complex-Gaussian RV with zero mean and unit variance. 
The normalized correlation coefficient between h p ji (r + r') and hji (r + r') satisfies p Pj . = \J~vp\J~&^.- 

For simplicity, the notation h S) ji is used to represent hji (r) and h Pi ji (r + r'). Specifically, for the RS scheme 
((5]), h s ,ji = hji(r), and for the RS scheme ©, h s> ji = h Pi ji(r + r'), then 7 SJ j = |/i Si jj| 2 . With the unified 
notation /i s ji, the analysis of the RS schemes §5§ and © can be expressed in a unified manner. 



IV. Ergodic Capacity of Analysis 

To analyze the ergodic capacity of bidirectional AF RS in ©, the distribution functions of % jk = hjk ( T + T ' 
in © need to be obtained first. 

Lemma 1: According to the RS schemes and the relationship between h t ji and h s ji, the PDF of ^t,jk is 
expressed as 



N N-l 



A,* (*) = ££ £(-!)* i + E 



i=l t=0 A t 





cxp 




where 




4 . *V 



+- 



.7' 



-i 



.7' 



2 2 / . 2 




(7) 



(8) 



(9) 



In addition, Y^a * s tne abbreviation of A t c{i,...,N}\i > an d l-^tl represents the cardinality of set A t . Moreover, 



\A t \=t 



Ug ji and afj i are the variances of h s ji and ht,ji, respectively. Specifically, for the RS scheme ©, o\ ^ = cr? , 
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T t,Jt 



<i- ., a? , pji = peP^, and p e = o 2 ^ Jo\.. Q; for the RS scheme ©, a 2 sji = a 2 . t , a 2 ^ = cr| , and 



Pii = PeP Pji - AISO, O-^ = ^ ^^^./(o-^. + a 2 2 J_ 

Proof: The derivation is given in Appendix B. 
Lemma 2: 



(a, m, n) 



In (y + a) e 



-m(y— x) — nx 



dxdy 



tp(a,n)-<p(a,m) , 

m-n ' T <h 

l—ame am E 1 (am)+mip(a,m) 



(10) 







,m = n. 



where ip (a, 6) = / In (x + a) e^dx = ^a+e^iab) ^ ^ 331]i and El (x) is the exponential integral function 

o 

na. 

Proof: Lemma 2 can be achieved by applying the integration by part, then discussing under the situations 
that m = n and m ^ n. ■ 
Proposition 1 : Applying the Lemma 1 and Lemma 2, the ergodic capacity of bidirectional AF RS is 



^ _ t\ + r 2 - r 3 - r 4 

21n2 



(11) 



where 



AT AT-1 



<7 



j=l i=o A t 



(-1)* 



/GA t <M 



0V V A^t,ji 



t,ji 
1 



—99 m, - — T 

Ipr V f ff t 



+—99 n 

V>s \ ^til 



N N-l 



N N-l 



<-j+2 



EEEEEE 



-i)'(-i 



■0 S + lp r J -0s 4=1 t=0 A t i'=l t'=0 A t , °~' J* ^ J» 
8 I ran, — '4-, — I +G© ran, — 




(12) 



1+E 



a 



2 

s,ji' 



at 



Ips+tpr 0s y 



+0Cj© ( mn, 



0? s +0 r 0s / 



(13) 



s + VV 0s y \ 0s + 0r 0s J 

In addition, m = ift s /ijj s , and n = (*0 S + tp r )/tp r . and £j can be obtained by £ 3 - and Q, respectively, by 
substituting i,l,A t in ([D and Q with i',l',A t >, respectively, and {j,j} = {1,2} or {2, 1}. 

It is noted that for the symmetric network structure, i.e., a 2 = 1, pji = p, i = 1, . . . ,N, j = 1, 2, we have 

= ( N ^ 1 ), and ^TJ^/ = (^T 1 ), thus the expression of capacity in Proposition 1 can be further simplified. 
Proof: The derivation is given in Appendix C. ■ 



V. Simulation Results and Discussion 

In this section, Monte-Carlo simulations are provided to validate the preceding analysis and to highlight the 
performance of bidirectional AF RS with outdated CSI. Without loss of generality, the sources and the relays 
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are assumed to have the same transmit powers, i.e., p s = p r = Pq. The network structure is assumed to be 
symmetric, i.e., a\_. = 1 and f dji = f d , j = 1, 2, i = 1, . . . , N. 

Fig. 1 studies the impact of outdated CSI on the ergodic capacity when N = 4 and adopting the RS scheme 
(f5]). The x-axis is SNR = Pq/u 2 , and the imperfect estimation is considered, i.e., the variance of detection error 
a e = v'n/Po ifTOl . Different lines are provided under different f d T, where larger f d T means CSI is outdated 
more severely than smaller f d T, and f d T = means CSI is not outdated. The figure verifies that the expression 
of Proposition 1 is the tight lower bound of the simulated results. We also observe that the ergodic capacity 
degrades once CSI is outdated, e.g., the performance loss between fdT = 0.3 and f d T = is about 4 dB in high 
SNR, and more severely outdated CSI results in greater performance loss, although outdated CSI has no impact 
on the multiplexing gain. 

Figs. 2-4 study the capacity of bidirectional AF RS with outdated CSI, when adopting the RS scheme ((6]) 
and the variance of detection error a 2 e = 0. We assume the time interval between relay selection process and 
its subsequent data transmission process satisfies t' = 1, and the time interval between adjacent relay selection 
processes satisfies A = 2. The prediction lengths of different channels are assumed to be the same, i.e., Lj% — L, 
J = 1,2, i = l,...,N. 

Fig. 2 plots the simulated and the analytical ergodic capacity of outdated CSI versus SNR = P /al, when the 
coefficient of outdated CSI is fixed, i.e., f d T = 0.3. The simulated results verify the lower bound in Proposition 1 
is tight when the channel prediction is adopted. For the symmetric network, the line of L = 1 represents the 
RS scheme (f5]), because p p = pf when L = 1. Also, the line of L = oo means the CSI is perfect, because 
p p = 1 when L = oo. As this figure reveals, the RS scheme based on channel prediction © outperforms the 
scheme without channel prediction d5j, and the performance gain gets larger when increasing the prediction 
length. The qualitative explanation of the phenomenon is that increasing the prediction length results in growing 
the correlation coefficient p p , thus the prediction becomes more accurate and the performance gets improved. 

Fig. 3 investigates the impact of f d T on capacity when SNR = 10, 15, and 20 dB. As the figure reveals, the 
curves of capacity under different prediction lengths L have almost the same performance in small f d T, and all 
the capacity degrades as f d T increases. However, larger L has better robustness of capacity, e.g., for L = 1 and 
SNR = 15 dB, capacity degrades to 3 bps/Hz when f d T = 0.18, whereas for L = 2, it is f d T = 0.27 when 
capacity degrades to the same level. Therefore, larger L improves the robustness of capacity. 

Fig. 4 investigates the impact of prediction length L on the capacity. The y-axis is the normalized difference 
of capacity, i.e., the difference between the capacity obtained by the RS scheme (O and the capacity obtained 
by the optimal RS scheme (@]), normalized by the latter. Smaller normalized difference means the capacity of the 
RS scheme (O with outdated CSI has closer performance with the optimal performance. Different lines of Fig. 4 
are plotted under different f d T and SNR, and all the normalized differences decrease monotonously to zero as L 
increases, which reveals that the RS scheme © can infinitely approach to the optimal performance, by enlarging 
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the prediction length. Furthermore, the results also reveal that the lines with larger f^T needs larger L to satisfy 
the same need of the normalized difference. 

VI. Conclusions 

The impact of imperfect CSI on the ergodic capacity of bidirectional AF RS has been investigated in this 
paper. The tight lower bound of ergodic capacity is derived in a closed form and verified by simulations. The 
results reveal that imperfect CSI will jeopardize the capacity of bidirectional RS network, whereas the RS based 
on the channel prediction can compensate the performance loss, and it can infinitely approach to the optimal 
performance, by enlarging the prediction length. 

Appendix A: Explanation of the Optimal Scheme 

To maximize the ergodic capacity of © is equivalent to maximize 71*72*, according to the approximation 
log 2 (l + 7ji) ~ log 2 7**, in which 7^ is provided in ©. Using the approximation xy/(x + y + c) « min (x, y), 
7ii72i is bounded by min (|/^,ii| 2 , \ht,2i\ 2 ) 2 - Therefore, (J3j is optimal in maximizing the capacity. 

Appendix B: Proof of Lemma 1 



Similar to [6], the PDF of jt,Xk can be expressed as 

d N 

Am* ( z ) = ~IZ Pr < zr]k = i} 



(b) 



8=1 

N 00 



Y / A*,i«i7«,» A.,» (y) Pr {7*,ii ^ %M%,u = y} Pr i k = *I*t»,x* = y, %,u < dy 

1 h 

00 

+ 5^/Am<I^,h (^ly) A.,i* (y) Pr {7s,ii>7s,2i|7s,ii = y}Pr{fc=«|7 Si i i =y,7 Si ii>7 Sj 2i}dy (14) 

where (a) is satisfied by the total probability theorem, which divides the union event jt,ik < z into iV disjoint 
events, i.e., E4 is the best relay and 7^1* < z, i = 1, . . . , N; (b) is fulfilled by the division of two disjoint events, 
i.e., 7 S) ij > 7 S 2t and 7 Sj ij < % 2i- Furthermore, I\ and I 2 can be obtained by the definition of the RS schemes 
((5]) and ((6]), the order statistics, and lITTl eq. 26]. Substituting the exponential distributions of 7 S) i* and 7 Sj 2i into 
(fT4l . Lemma 1 is verified by applying the conditional PDF of it,lk an d 7s,lfe and lfl4l eq. 6.614.3]. Moreover, 
the PDF of 7 t 2 fc can be obtained similarly. 
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Appendix C: Proof of Proposition 1 
The ergodic capacity of bidirectional RS is tightly bounded by 



_(«) 
C > E. 



1 



'7t,ifc,7t,2fc 



log 2 



Wv7t,i*+ m ) {ip s %2k+ n) 



+ E. 



7t,ife)7t,2fc 



log 2 



W>s7tifc+ n ) (A%2k+ rn) 



(15) 



*/>s7t,i*+ ( fps+ipr ) lt,2k+mn 
where m = ipg/ipg, and n = (ij) s + ip r ) /ip r . Comparing with d2), (a) is achieved by adding the constant mn — c 
in the dominator of (fT5l . which has little effect on the performance, like the SER analysis 0. 

After some manipulation of the logarithmic and expectation operations, the capacity of (fl3T ) is rewritten as 
C = Tl+T 2\n2~ Ti ' wnere T\ = E^ 4 lfc {in {^rlt,ik + w)+ln (V ; s7t,ifc+ n ) }> and T 2 can be obtained by substituting 
%ik in T\ with 712ft- According to the fact that fy (z) = fx (z/m) jra when Y = mX (m > 0), the PDFs of 
ipsjtjk an d ipr'Jtjk can be obtained by the PDF of j = 1,2, then (fT2l ) can be achieved by Lemmas 1 and 



2. Moreover, T3 



Jt,lk ,7t,2fe 



In 



, and T4 can be obtained by permuting ^t,ik 



lp a +1pr J A ft,lk+1ps A ft,2k+ mn 
with j t ,2k in T3. j t ,ik and 7 ti 2fc are independent RVs, thus the PDF of {^ip s + ip r j jt^k + ■0 s 7t,2fc can be achieved 
by the convolution of n/S s + ?/) r J gift's and Tpsjt^k's PDFs, then T3 in (fT3T ) can also be achieved by Lemmas 1 
and 2. The expression of T4 can be obtained similarly. 
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Fig. 1. Capacity under outdated CSI when applying the scheme Q. 




Fig. 2. Capacity when applying the RS scheme {6}. 
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Fig. 4. The normalized difference of capacity versus prediction length L. 
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